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ABSTRACT 

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame 
are examined mathematically by a direct analysis of the Navier-Stokes equations. It is 
proven that two such shear flows are dynamically similar if and only if their initial di- 
mensionless energy spectrum E*(k*, 0), initial dimensionless shear rate SK 0 /so, initial 
Reynolds number K^/veo, and the ratio of the rotation rate to the shear rate ti/S are 
identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, 
they will only depend on the single parameter f 1/5. The commonly assumed dependence of 
such equilibrium states on fl/S through the Richardson number Ri — — 2(0/5')(l — 20/5) 
is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more 
than a weak approximation. To be more specific, Richardson number similarity is shown to 
only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., 
to certain second-order closure models) wherein closure is achieved at the second-moment 
level by assuming that the higher-order moments are a small perturbation of their isotropic 
states. The physical dependence of rotating turbulent shear flows on fl/S is discussed in 
detail along with the implications for turbulence modeling. 
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1. INTRODUCTION 


A variety of geophysical and engineering systems of interest contain turbulent shear 
flows which evolve within a rotating framework (e.g., ocean currents, atmospheric boundary 
layers, turbomachinery, etc.). Consequently, it is essential that the effects of rigid body 
rotations on turbulent shear flows be understood if such physical systems are to be modeled 
properly. In order to gain insight into these complicated systems, a variety of researchers 
have considered the idealized problem of plane homogeneous turbulent shear flow in a 
steadily rotating frame (see Bradshaw 1 , Ferziger and Shaanan 2 and Bardina, Ferziger, and 
Reynolds 3 ). This work, which is somewhat empirical in its orientation since it is based 
on turbulence modeling, has shed some important new light on the physical structure 
of rotating turbulent shear flows. However, several fundamental questions related to the 
scaling of rotating turbulent shear flows remain unanswered. To be specific: (a) which 
dimensionless parameters must be invariant to ensure the similitude of two homogeneous 
turbulent shear flows in a rotating frame, (b) when will equilibrium states exist that 
depend only on the ratio of the rotation rate to the shear rate f J/5, and (c) when is this 
dependence on f 1/5 exclusively through the Richardson number Ri — (2fi/5) (1 — 20/5)? 
The answers to such scaling questions play a crucial role in the construction of turbulence 
models that have the correct physical behavior in rotating frames. The motivation for the 
present paper is therefore to derive scaling laws for homogeneous turbulent shear flows in 
a rotating frame based on a direct analysis of the Navier-Stokes equations. 

It will be shown that two homogeneous turbulent shear flows in a rotating frame are 
dynamically similar if and only if their initial dimensionless energy spectrum, shear rate, 
and Reynolds number (all scaled by the initial turbulent kinetic energy and dissipation 
rate) are identical in addition to the ratio of the rotation rate to the shear rate Cl/S. 
It then follows that if an equilibrium structure exists at high Reynolds numbers which 
attracts all initial conditions (i.e., a universal equilibrium), its turbulence statistics will 
be determined by a single parameter - the ratio of the rotation rate to the shear rate 
Cl/ S. We will examine the circumstances under which the dependence of such equilibrium 
states on Cl/S collapses to Richardson number similarity. It will be proven that Richardson 
number similarity does not hold for general solutions of the Navier-Stokes equations for 
homogeneous turbulent shear flow in a rotating frame. However, it will be demonstrated 
that low order truncations of the Navier-Stokes equations such as certain second-order 
closure models (which are obtained by assuming that third-order moments are a small 
perturbation of their isotropic state) do exhibit Richardson number similarity with respect 
to the turbulent kinetic energy, dissipation rate and Reynolds shear stress. This result 
has important implications for the modeling of rotating turbulent shear flows where it has 
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been a common practice to modify simpler turbulence models by the inclusion of an ad 
hoc Richardson number dependence to account for rotations 4-6 . The broader implications 
of these results for the modeling of rotating turbulent shear flows will be discussed briefly 
in the last section. 
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2. DIMENSIONAL INVARIANCE OF ROTATING TURBULENT SHEAR 
FLOWS 


We will consider the problem of homogeneous turbulent shear flow in a rotating frame 
for a viscous incompressible fluid (see Figure l). The governing equations of motion for 
this flow are the Navier-Stokes and continuity equations given by 7 

dv 

— — b v • Vv = — VP + r/V 2 v — 2H x v (l) 

at 

V • v = 0 (2) 

where v is the velocity field, P is the modified pressure (which includes the centrifugal 
acceleration potential), v is the kinematic viscosity of the fluid, and 

n = nk (3) 

is the rotation rate of the reference frame relative to an inertial framing. As in the usual 
treatments of turbulence, the velocity field v and pressure P are decomposed into ensemble 
mean and fluctuating parts as follows 

v = v + u, P = P + p. (4) 


For those turbulence statistics which are homogeneous, the ergodic hypothesis can be 
invoked and ensemble averages can be replaced by spatial averages. The spatial average 
of a flow variable / is defined as follows 8 


f^ = ^vlv f{x ’ t)dl 

For the plane turbulent shear flow under consideration, 

dvi 
dxj 

or, equivalently, in Cartesian tensor notation 

dvi 


x. 


( 5 ) 


0 

s 

0 ^ 

0 

0 

0 

0 

0 

0 J 


( 6 ) 


dxj 


ss u s 2i 


( 7 ) 


where S is the shear rate and is the Kronecker delta. It is a simple matter to show that 
the mean velocity field given by IT,- = SSuS^jXj is consistent with the Reynolds equation 
and mean continuity equation 


b v • Vv = — VP + z/V 2 v -20 xv 

dt 


( 8 ) 


3 



V • v = 0 (9) 

which are obtained by taking the ensemble mean of Eqs. (l)-(2) and making use of the fact 
that the divergence of the Reynolds stress tensor r,-y = —ujuj vanishes since the turbulence 
is homogeneous. Specifically, v identically satisfies the continuity equation (9) and the 
Reynolds equation (8) if 

P = — 0 Sy 2 + constant. (10) 

The fluctuating velocity u is a solution of the momentum and continuity equations 
given by 

dn 

— + v • Vu = — u • Vu — u • Vv — Vp + uV 2 u -20xu (11) 

at 

V • u = 0 (12) 

which are obtained by introducing the decomposition (4) into (l)-(2) and then differencing 
the resulting equations with (8)-(9). For the rotating homogeneous turbulent shear flow 
under consideration, equation (ll) takes the component form 


dUi . o* * du i 


du dx) 

~ Uk ~faT ~ ~ (^*'1^*2 — 2S7er 3t - k )uk + uV 2 Ui (13) 


where is the permutation tensor. 

Now, we will non-dimensionalize these equations by introducing the dimensionless vari- 
ables 

* u i * x i £ 0 /, 


f = st (15) 

where K 0 and e 0 are, respectively, the initial turbulent kinetic energy and dissipation rate. 
It should be noted that since the turbulence is homogeneous, the turbulent kinetic energy 
K and dissipation rate e defined by 8 

. , 1 . , dUi du; , , 

KW = 2 UiU " s{t)= v s7 i d7 i (16) 

are functions of time alone. The introduction of (14)-(15) into (12)-(13) yields the dimen- 
sionless equations of motion 


du i , c c * du i 

~p>+* + 


, du* { dp* 
SK 0 Uk dxi dx* 


A * * 

Aik u k 


SKq Re o 


V*V 
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where 


(19) 


Re 0 = — 2- 

U£q 

0 1 - 2fl/5 0 \ 

A; k = 6 a S k2 -2{Q/S)e 3ik = 20/5 0 0 

l 0 0 0 / 


The equations of motion (17)-(18) must be solved subject to the initial condition 

u*(x*,0) = u^ 0) (x*) (21) 

at time t* = 0. We will restrict our attention to the case where the turbulence is initially 
isotropic - the same initial condition that has been used in previous physical and numerical 
experiments on homogeneous turbulent shear flows (c.f., Champagne, Harris, Corrsin 9 , 
Tavoularis and Corrsin 10 , Rogallo 11 , and Bardina, Ferziger, and Reynolds 3 ). Given an 
isotropic turbulence, the hierarchy of turbulence statistics are derivable from the energy 
spectrum E(k,t). For example, the turbulent kinetic energy K and dissipation rate e are 
given by 

roo 

K{t) = / E(k,t)dk (22) 

Jo 

e[t ) = 2 u [ k 2 E(k,t)dk (23) 

Jo 

and the skewness S k in the quasi-equilibrium range 12 takes the form 

* = L k ^L (24) 

7 [J k 2 E(k)dk]* 

Hence, if we are only interested in the time evolution of the Reynolds stress tensor r t y and 
turbulent dissipation rate e, it is reasonable to believe that we need only be concerned with 
initial conditions on the fluctuating velocity u*(x*,0) through the initial energy spectrum 
E*(k*,0). These initial fields are related by 


1 roo 

E*(k*,0) = — F 2 J (5,*,-(r*)e -4k *' r *d 3 


where 


Q*ij( r*) = u?(x*,0)u*(x* + r*,0) (26) 

3 

is the dimensionless two-point velocity correlation tensor and k* = kK§ fe is the dimen- 
sionless wavenumber. 

Thus, it has been demonstrated that two homogeneous turbulent shear flows in a 
rotating frame (evolving from an initially isotropic state) are dynamically similar in a 
statistical sense if the following flow parameters are identical: 
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• the initial dimensionless shear rate SK 0 /e 0 

• the initial dimensionless energy spectrum E* ( k * , 0) 

• the initial Reynolds number K\Jveq 

• the ratio of the rotation rate to the shear rate fi / S. 

When these four parameters are the same for two different homogeneous turbulent shear 
flows in rotating frames, their dimensionless Reynolds stress tensor r/y and dissipation rate 
e* (as well as other moments constructed from u*) will evolve identically with respect to 
the dimensionless time t* = St. Consequently, if universal equilibrium states exist (which 
are solutions for t — ► oo that are completely independent of initial conditions) , they will 
only depend on a single parameter - the ratio of the rotation rate to the shear rate U/S. 
In the next section, we will explore when such universal equilibrium states depend on Cl/S 
through the Richardson number — 2(0/S)(l — 2f 1/S). 
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3. RICHARDSON NUMBER SIMILARITY 


Bradshaw 1 presented an analogy between rotating and stratified turbulent flows. In 
that study, he defined a Richardson number for rotating shear flows as follows 

-20(5-20) 


Ri = 


5 2 


(27) 


This is analogous to the Richardson number for density stratified flows in that (27) serves 
as the square of the Brunt- Vaisala frequency ubv for small amplitude inertial oscillations 
induced by the application of a rotation to shear flow. Although the analysis leading to 
Eq. (27) is straightforward, we feel that it would be instructive to include it here since 
Bradshaw 1 did not provide the details. If we consider small homogeneous perturbations of 
shear flow of the form 

v = 5yi + u(£) (28) 

P = - nsy 2 + p(t) (29) 

it is a simple matter to show, by substituting Eqs. (28)-(29) into (13), that u a solution 
of the equation 

ilk = —AktUi (30) 

where A kt = S8 kl 8 t2 — 2tte Ske . Equation (30) has a general dimensionless solution of the 
form 


«; = a 1 e‘^ m '• + r 

(31) 

u, = a' 1 e , ' /m ■’ + *' 

(32) 


where a,- and a\ are related to the initial conditions. Hence, the Brunt-Vaisala frequency 
for this perturbation of homogeneous shear flow in a rotating frame is given by 


U>BV — VRl 


(33) 


as argued by Bradshaw. This motivated Bradshaw 1 to propose a modification of the mixing 
length, to account for rotations in turbulent shear flows, that is of the form 

^ = i - m. (34) 

In (34), l is the mixing length with rotations, io is the mixing length in the absence of 
rotations, and j3 is an empirical constant. Various investigators in the intervening years 
have modeled rotating turbulent flows with a mixing length of the same form as (34) (c.f., 
Lezius and Johnston 4 , Launder, et al. 5 , and Howard, et al. 6 ). When such a length scale is 
used in a turbulent eddy viscosity model, the Reynolds shear stress, as well as the turbulent 
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kinetic energy and dissipation rate, will exhibit Richardson number similarity. Hence, it is 
critical to know in what sense such similitude is consistent with the Navier-Stokes equations 
if we are to properly evaluate these turbulence models. 

We will now demonstrate that Richardson number similarity is not consistent with 
general solutions of the Navier-Stokes equations. In order to accomplish this task, we will 
eliminate the pressure from equation (17) by taking the divergence of that equation and 
making use of the constraint (18). This yields the Poisson equation 


e 0 du* k du ? 
SK o dxj dx* k 


which has the solution 


{A*ik + £*' 1 ^* 2 ) 



(35) 


V 


/: 


£ 0 du* k du\ 
An r J-oo |x* - z*\ [SKo dz ? dz k 

in an unbounded flow domain where 

B*ik ~ A ik + f>nf>k 2 


, R* du * k 


3 * 


d*z 


(36) 


0 2-2 n/s 0\ 

2 n/s o o 

^ 0 0 0 J 


(37) 


Hence, u* is a solution of the integro-differential equation 


du* 


du*. 


dr +SllS * x ’idx 


* du* 

U u 


^0 


SK 0 k dx* k 


A >‘ + SK„Re r 


7*2 




1 d 


47 r dx 


■/: 


^0 


du* k du\ 


S K 0 dz*, dz 


B 


Ik 


du l 

dz*. 


3 * 


d z 


(38) 


Since the eigenvalues of A*^ are ±i\/~Ri and the eigenvalues of B*- are Bi\J Ri — 20 / S 
there exists no basis relative to which the dependence of (38) on f l/S can be collapsed to a 
Richardson number dependence. Hence, Richardson number similarity does not follow from 
the full Navier-Stokes equations. The recent large-eddy simulations of Bardina, Ferziger, 
and Reynolds 3 were the first to cast grave suspicions on such Richardson number scaling. 
In figure 2, results of these large-eddy simulations for homogeneous turbulent shear flow 
in a rotating frame are shown for 11/5 = 0 and f l/S = 0.5 - the two values of fi/5 
which correspond to a Richardson number Ri = 0. Clearly, the turbulent kinetic energy 
and dissipation rate associated with these two cases are quite different and thus indicative 
of a violation of Richardson number similarity (the initial conditions for both flows were 
identical) . 
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We will now proceed to show that certain low-order truncations of the Navier-Stokes 
equations do exhibit Richardson number similarity. Second-order closure models will be 
considered wherein the infinite hierarchy of moments obtained from the Navier-Stokes 
equations are truncated at the second-moment level by modeling the higher-order moments. 
This second moment level consists of the Reynolds stress transport equation 


i~ij — -Ai k^kj -AjicTki 


(39) 


where r t y = —UiUj is the Reynolds stress tensor and given by 


= -p 



+ 



+ 2 v 


dui duj 
dxk dxk 


(40) 


is the sum of the pressure strain and dissipation rate correlations. Equation (39) (which is 
not a closed system for the determination of r tJ ) is obtained by taking the ensemble mean 
of the symmetric part of the product of equation (13) with Uj. Closure can be achieved 
by assuming that 3>,y is a functional of the variables r, Vv and e as follows in an inertial 
frame* 

$ij = $«/(b, Vv, K, e) (41) 

where = — (r,y + | K6ij)/K is the anisotropy tensor. Since for isotropic turbulence, $>,-,• 
is given by 8 

$«'; = \ eS H* ( 42 ) 

the physical motivation for (41) is clear: anisotropies in arise from the mean velocity 
gradients Vv which are then reflected in nonzero values of b. For small perturbations of 
isotropy, we have 

b = (5 b, Vv = <5(Vv) 


To the first order in <5b and <$(Vv), is given by 

'dSu' 


<Mb,Vv,K,e) = ($,y)o + 


■ *j 


dh 


(5b + 


V^Vvj 0 


<5(Vv) 


(43) 


2 , _ _ / dv{ dvj 

+ c.eki + c 2 K + — 


where (-) 0 denotes a function evaluated at b = 0, Vv = 0 (i.e., for the state of isotropy) 
and Ci and c<i are dimensionless constants. Equation (43) is obtained by making a Taylor 
expansion about the isotropic state (see the Appendix). The modeled equation for 3> tJ - is 

*The reader should note that in a rotating frame dvi/dxj must be replaced by dvi/dxj + (see 

Speziale 12 ). 
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precisely the same as that used in the Rotta-Kolmogorov model 13 if we take ci = 3.21 and 
c 2 = -0.112. 

As an alternative to this approach, the pressure strain and dissipation rate correlations 
can be modeled separately. Since dissipation is manifested primarily at the small scales 
(which are not far removed from isotropy at high Reynolds numbers) the dissipation rate 
correlation can be approximated by its isotropic form 


2v- 


tj- 


(44) 


du{ duj _ 2 
dx k dxk 3 

The pressure strain correlation vanishes for isotropic turbulence and for homogeneous 
turbulence can be shown to be of the form (see Shih and Lumley 14 ) 

dv k 


( dui duj \ 

p (aij + aiy = + Mii " 


dxt 


(45) 


in an inertial framing where R%j{t) and M{jki(t) are functionals of the energy spectrum 
tensor Eij(k,t) over wave number space. It can be assumed that (see Shih and Lumley 14 ) 


Rij = Rij(b,K,e) (46) 

M{j ke = Mij ke (b, K, e ) (47) 


For small anisotropies, (46)-(47) can be expanded in a Taylor series which, to the first 
order in b, yields the following expression for the pressure strain correlation in a rotating 
frame 15,16 : 


~P 


du{ duj \ 
dxj ^ dxi J 


+ C 2 K 


dvj dvA 
dxj dx, ) 


+ CzK(bi k S j k + bj k Si k 


-^buSuSij) + C 4 K{b ik W jk + b jk W ik ) 


(48) 


where 


*' 2 \dx, Ox, I ' ” 2 [ dXj Ox,) + ” 

and we have made use of the fact that the pressure strain correlation is traceless and must 
vanish in the limit as b, Vv, and Cl go to zero. If we take C\ = 1.8, C 2 = — 0.4, C 3 = —0.6, 
and C 4 = —0.6, the Launder, Reece, and Rodi 17 model is obtained. 

When closure relations for of the general form (41) are substituted into (39), closure 
is achieved once a modeled transport equation for e is provided. For high Reynolds number 
turbulence that is homogeneous, the exact transport equation for £ takes the form 17 



De dui dui du k , d 2 u,- d 2 u,- 

= —2v 2u 

Dt dx k dx t dxi dx k dx t dx k dxi 


(50) 
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( 51 ) 


It can be shown that the lowest-order model for (50) is of the form 

De e e 2 

Di = Cel K P ~ Ce2 K 


where P = Tijdvi/dxj is the turbulence production and C e i and C e 2 are dimensionless 
constants (see Lumley 18 and Launder, Reece, and Rodi 17 ). Equation (51) is obtained by 
assuming that the higher-order correlations on the right-hand-side of (50) are functions of 
b, Vv, K, and e which are expanded in a Taylor series to the first order in b (e and K merely 
set the length and time scales for dimensional consistency). For the Rotta-Kolmogorov 
model, 5 C e i = 1*8 and C e i = 2.0 whereas for the Launder, Reece, and Rodi Model, 
C e i = 1-45 and C e2 = 1.90. When (51) is solved in conjunction with (39) and a model 
of the form of (41), a closed system of equations for the determination of the Reynolds 
stress tensor and the turbulent dissipation rate are obtained. This constitutes a low order 
truncation of the hierarchy of moments of the Navier-Stokes equations (achieved at the 
second moment level) by assuming that higher-order moments are a small perturbation of 
their isotropic state. 

Equation (39) can be written in the alternative dimensionless matrix form 


7 * = — A* 7 * - 7 * A* 


(52) 


where r* = t/K 0 , $* = $/K 0 S and A* T denotes the transpose of A*. If we introduce the 
variable 

t** = A V + 7 * A* t , (53) 

then (52) can be written in the alternative form 

r * = r* - av* - 7 **a* t (54) 


where 

$** = A*$* + $*A* t (55) 

Equation (54) is obtained by adding the two equations that result from the pre-multiplication 
of (52) with A* and the post-multiplication of (52) with A* T . Since, 

0 1-2 n/s 


A* = 


it follows that 


2U/S 

2(1 - 2n/S)r 1 * 2 


4(n/sK 2 J 


(56) 


(57) 


§The Rotta-Kolmogorov model is actually based on a transport equation for the length scale which, for 
homogeneous turbulence, can be converted to a transport equation for e of the form (51). 
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AV 


* * a *T 

“= T A — 


(1-20 /5)r* -2Ri r* 


—2Ri r- 


12 


20/5 


12 


'd J 


where 


= (1 - 2CI/S)t22 + 20/5 r* x 

Taking the trace of (54) and dividing by two we obtain the equation 

1 


*« = > r * - 


whereas the 1,2 - component of (54) yields the equation 


*2 = $ 12 - 2 Ri T* 2 . 


(58) 


(59) 


( 60 ) 


(61) 


Equations (60) -(61) must be solved in conjunction with the turbulent kinetic energy and 
dissipation rate equations given in the dimensionless form 


K* = r* — e ° 

12 SK n 


*2 


. * £ * s 0 e 

' ~ C ‘ 1 K :Tn ~ C ‘ 2 SK 0 K 7 


(62) 

(63) 


where if* = K/Kq and e* = e/eo- The system (60)-(63) is solved subject to the initial 
conditions 

JT(0) = 1, e*(0) = 1 (64) 

2 


ri,(0)=0, rj(0) = 


(65) 


Hence, for a given initial value of eq/SKo, (60)-(63) will yield solutions for K*,e*,t1 2 and 
r* d which exhibit Richardson number similarity provided that 


tr$** = /(n*2, r hK*, s*,e 0 /SKq, Ri) 
$11 = g(rl 2 ,T*,K*,E\E 0 /SK 0 ,Ri ) 


( 66 ) 

(67) 


where / and g are arbitrary functions that need only be sufficiently smooth. For the 
Rotta-Kolmogorov model, 


I ** 2 ^0 * £ 0 £ TS* f *\ 

$»=--=* -C 177 ^ — {-K +T d ) 


3 S Kn 


' SK 0 K* 3 


( 68 ) 

(69) 
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which satisfies (66)-(67). Hence, the Rotta-Kolmogorov model exhibits Richardson number 
similarity. However, in the Launder, Reece, and Rodi model, 


tr $** = -2 C x 


So 


-t* 2 + 2 C 2 K* - 2C S K* 


+ 


SK 0 K* 

2C,(1 - 2fi/S) (|i- - ^ 


4 Ti 1 
- + — + 
.3 K* 


(70) 


*** 

*12 


2 £o 
3SKq 


£* - Ci 


so 


SK 0 K' 




2n\ 

T) T 


12 


(71) 


which clearly violates (66)-(67). Consequently, the Launder, Reece, and Rodi model vio- 
lates Richardson number similarity. In figure 3, computed results for the time evolution of 
the turbulent kinetic energy is shown for a Richardson number Ri = 0 (which corresponds 
to the two dimensionless rotation rates of 0/5 = 0 and 0/5 = 0.5) for the second-order 
closure models considered in this study. Consistent with the proofs already presented, 
the Rotta-Kolmogorov model yields the same results for both cases whereas the Laun- 
der, Reece, and Rodi model yields results that are dramatically different (and which, in 
comparison to the large-eddy simulations, deviate too strongly from Richardson number 
similarity). Both the Launder, Reece, and Rodi model and the Rotta-Kolmogorov model 
predict a universal equilibrium for the dimensionless shear rate SK/e and anisotropy ten- 
sor bij in the limit as f -> oo. These results are shown in Table 1. Clearly, the results 
obtained from the Rotta-Kolmogorov model exhibit Richardson number similarity whereas 
those obtained from the Launder, Reece, and Rodi model do not. 

Finally, we would like to comment in more detail on the asymptotic long time be- 
havior of the Rotta-Kolmogorov model which constitutes the lowest-order second moment 
truncation of the Navier-Stokes equations. It can be shown that for dimensionless time 
t* » 1, the Rotta-Kolmogorov model yields solutions such that (see Speziale and Mac 
Giolla Mhuiris 19 ) 


K*(t*) - exp 


(612)00 + 


SK 


(72) 


**(*•) 


exp 


(£> 12)00 + 


SK 


(73) 


for intermediate rotation rates —0.11 < Cl/S < 0.61. For values of 0/5 outside of this 
range, the flow undergoes a relaminarization (i.e., K and e — ► 0 as t — > 00 ). In figure 4, 
the time evolution of the turbulent kinetic energy computed from the Rotta-Kolmogorov 
model is shown for rotation rates in both characteristic regimes (i.e., for 0/5 = 0.25 in the 
unstable regime and for 0/5 = 0.75 where a relaminarization occurs). It is interesting to 
note that the existence of an unstable flow regime, at intermediate rotation rates, where 


13 



there is an exponential growth in time of the turbulent kinetic energy and dissipation rate 
has been argued based on alternative models (see Bertoglio 11 , Rogallo 20 , and Tavoularis 21 ). 
Furthermore, the Rotta-Kolmogorov model predicts the existence of universal equilibrium 
states in this unstable regime where the anisotropy tensor and dimensionless shear rate 
approach the values of (&»y)oo and ( SK/e )<» which depend only on Cl/S through the 
Richardson number 19 - a result which is a crude approximation of numerical simulations 
of the Navier-Stokes equations 3 . 
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4. CONCLUSION 


Conditions for the dynamic similarity of homogeneous turbulent shear flows in a rotat- 
ing frame have been obtained based on a direct analysis of the Navier-Stokes equations. 
It was proven that the dimensionless parameters which established similitude are the ini- 
tial energy spectrum !£*(£*, 0), the initial shear rate SKo/eo, the initial Reynolds number 
K^/veo, and the ratio of the rotation rate to the shear rate Cl/S. It was proven that the 
commonly assumed dependence of the Reynolds stress tensor and dissipation rate on Cl/S 
through the Richardson number Ri = — 2(Cl/S) (1 — 2 Cl/S) is not a rigorous consequence 
of the Navier-Stokes equations. In fact, only the lowest-order second-moment truncation 
of the Navier-Stokes equations (i.e., the Rotta-Kolmogorov model) exhibits Richardson 
number similarity with respect to the turbulent fields K*, and e*. The second mo- 

ment closure of Launder, Reece, and Rodi deviates quite strongly from Richardson number 
similarity. 

The introduction of an ad hoc Richardson number dependence in the length scale for 
a one equation turbulence model can only yield a crude approximation of the turbulence 
structure for rotating shear flows. However, the commonly used turbulence models tend to 
be so empirical in nature that the errors introduced in making the assumption of Richard- 
son number similarity may be comparatively unimportant. There is, of course, a need 
to develop turbulence models which can more accurately predict the turbulence structure 
in rotating turbulent flows. We believe that such improved models should be based on 
physically motivated modifications of the Reynolds stress and dissipation rate transport 
equations, which rationally account for rotational strains, rather than on the introduc- 
tion of an ad hoc Richardson number dependence. The development and testing of such 
improved turbulence models will be the subject of a future paper. 
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APPENDIX 


The Taylor expansion (43) can be written in the form 

2 Q-y 

(t>> Vv, K, e) = -eSij + c^kebu + (Al) 

o oxp 

where we have made use of representation (42) for (3> t -y) 0 and also of the fact that (5b = b 
while <5(Vv) = Vv. In this form of the expansion, Ciju and diju are defined by 


Cijke = 


d-ijkt — 


d*ij 

dbw 


d$ii 


, dv 


M/o 


(A2) 


(A3) 


where Vk,t = dvk/dxz. We remind the reader that (-)o denotes a function evaluated at 
b = 0, Vv = 0; therefore c^t and diju are functions of K and e alone. 

Dimensional and tensorial invariance along with the symmetry of <f>ij and bij yields the 
following form for Cijki'- 

c ijkz{K , e ) = Ci -{Sikfijt + bjkSa) + c^ebijdki (A4) 

where Ci and c[ one dimensionless constants. Since bu is traceless it follows that 

c i jktbkt — c i ebij (A5) 


The same arguments can be used to show that 


dijke{K,e ) — c 2 — (Sikbje + djkdu) + c' 2 K8ij6 k ^ 


and, hence, 


. dv k v 

dijktr . — = c 2 K 
OXp 


' dvi dvj 
dxj da;,-. 


(A6) 


(A7) 


where c 2 is a dimensionless constant. The final form of equation (43) then follows. 
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Rotta-Kolmogorov 

Model 

S3 

"co 

II 

0 

n/s = 0.5 

Launder, Reece 
and Rodi Model 

0 

II 

Co 

cT 

"to 

II 

p 

bi 

(SK/e) oo 

3.71 

3.71 

{SK/e)* 

5.42 

OO 

( 612)00 

-0.337 

-0.337 

( 6 12)00 

-0.369 

0 

(^ll)oo 

0.483 

-0.241 

(^ll)oo 

0.381 

-0.240 

( 622)00 

-0.241 

0.483 

( 622)00 

-0.190 

0.320 

{r d /K) CO 

-0.426 

-0.426 

{r d /K) oo 

-0.477 

-0.427 


Table 1 . Equilibrium values of the Rotta-Kolmogorov model and the Launder, Reece, and 
Rodi model for a Richardson number Ri = 0 . 


18 




































LIST OF FIGURES 


Figure 1. 
Figure 2. 

Figure 3. 

Figure 4. 


Homogeneous turbulent shear flow in a rotating frame. 

Time evolution of the turbulent kinetic energy for a Richardson number 
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Figure 3. Time evolution of the turbulent kinetic energy for a Richardson number 
Ri = 0 and initial condition e 0 /SKo = 0.496: o Large-eddy simulations 3 , 

- - Rotta-Kolmogorov model, — Launder, Reece, and Rodi model. 


22 




Figure 4. Time evolution of the turbulent kinetic energy obtained from the Rotta- 
Kolmogorov model for an initial condition e 0 /SKo = 1 (fi/S = 0.25 and 
fl/5 = 0.75). 
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